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This talk reviews the progress made in the determination of the light quark masses using lattice QCD and QCD sum 
rules. Based on preliminary calculations with three flavors of dynamical quarks, the lattice estimate is rUs = 75(15) MeV, 
a tantalizingly low value. On the other hand the leading estimates from scalar and pseudo-scalar sum rules are 99(16) 
and 100(12) MeV respectively. The r-decay sum rule estimates depend very sensitively on the value of |Vus|. The central 
values from different analyses lie in the range 115 — 120 MeV if unitarity of CKM matrix is imposed, and in the range 
100 — 105 MeV if the Particle Data Group values for |Kis| are used. I also give my reasons for why the lattice result is 
not yet in conflict with rigorous lower bounds from sum rule analyses. 



1 Introduction 

Quark masses are important fundamental parameters 
of the standard model. Our ability to extract them 
from first principle calculations of QCD will signal the 
onset of quantitative control over the non-perturbative 
aspects of QCD. In this talk I will summarize the cur- 
rent status of the extraction of light quark masses. 
All results for quark masses will be in the MS scheme 
at scale 2 GeV. Most of the time will be devoted to 
a review of the lattice results. A brief summary of 
sum-rule calculations will also be presented. A more 
detailed version of this review is being prepared in col- 
laboration with Tanmoy Bhattacharya and Kim Malt- 
man 

Of the light quark masses, I will concentrate entirely 
on the strange quark mass. The reason is that, at 
present, estimates of the ratios from chiral perturba- 
tion theory 



2rn, 



•mu + Tfid 
m,. 
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0.553(43) 



fact that current lattice calculations indicate that the 
the whole interval m„ ^ — TOc can be handled by the 
techniques used for light quarks. Consequently, charm 
quark can be simulated on the lattice with small dis- 
cretization errors. Following this approach, the AL- 
PHA collaboration has presented, in the quenched ap- 
proximation, the ratio [|2j 



^ = 12.0(5) 



(2) 



The advantage of their calculation, based on the 
Schrodinger functional approach and using the fully 
0(a) improved theory, is that they convert lattice 
masses to renormalization group invariant masses us- 
ing a non-perturbative method, thus avoiding the 
problem of the perturbative determination of the 
scale dependent renormalization constant Z')^ at scales 
1/a = 2 — 4 GeV and the subsequent matching of 
the lattice to a continuum regularization scheme at 
these scales. They find m^(2GeV) = 1.301(34) 
GeV, which compares well with the recent estimate 
(2GeV) = 1.26(4)(12) GeV by the SPQcdR col- 



(1) laboration [E]. 



are more accurate than the lattice results. Using them 
to extract to^ and md, once ms is known, avoids some 
of the uncertainties due to chiral extrapolations in lat- 
tice calculations, and of ignoring electromagnetic ef- 
fects in the simulations. 

There is a recent result on the charm quark mass that 
deserves mention. I mention this to also highlight the 

* It is a pleasure to thank Patricia Ball, Jonathan Flynn, and 
Laurent Lellouch for the opportunity to present this review, and 
to Tanmoy Bhattacharya and Kim Maltman for many discus- 
sions. I also thank S. Aoki, T. Kaneko, C. Bernard, V. Lubicz, 
H. Wittig, and G. Schierholz for discussions of their calcula- 
tions. The work of R.G. was, in part, supported by DOE grant 
KA-04-01010-E161. 



2 Lattice QCD 

A self-consistent determination of quark masses us- 
ing lattice QCD will be equivalent to the validation of 
QCD as the correct theory of strong interactions. In 
ideal lattice QCD simulations we need to dial six input 
parameters in the generation of background gauge con- 
figurations and then calculate quark propagators on 
them. These are the gauge coupling g, and the masses 
for up, down, strange, charm, and bottom quarks. The 
top quark is neglected because it is too heavy and too 
short lived. If we had petaflop scale computers we 
could carry out simuations with three light and two 
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heavy flavors of dynamical quarks with masses tuned 
to roughly their physical values m„, m^, rjis, nic, and 
mi,. The inclusion of heavier quarks, charm and bot- 
tom, in the generation of background gauge config- 
urations, however, are unlikely to significantly affect 
the vacuum structure and can therefore be safely ne- 
glected in the update for much less than 9 (or 100) 
GeV^. To study their interactions, charm and bot- 
tom quarks are incorporated at the stage of calculating 
quark propagators on these background gauge config- 
urations. Thus simulations with three light dynamical 
flavors are the goal of lattice calculations. 

Correlation functions in Euclidean space-time, from 
which various properties of hadrons are extracted, are 
constructed by tying together these quark propaga- 
tors and gauge link variables in appropriate combina- 
tions. For example the hadronic spectrum and asso- 
ciated decay constants are extracted from two-point 
correlation functions with the appropriate quantum 
numbers. The masses of hadrons, at least of the sta- 
ble ones and of those with very narrow widths, are 
determined from the rate of fall-off of these correla- 
tion functions at large Euclidean time. If QCD is the 
correct theory, these should agree with experiments 
up to electromagnetic corrections which are neglected 
in current simulations. 

Unfortunately, today's computers are not powerful 
enough to carry out calculations with physical val- 
ues of up and down quark masses. Instead, quarks 
heavier than their physical values, i.e. in the range 
TTis/S — TTis, have been studied. Also, isospin breaking 
and electromagnetic effects have been neglected, i.e., 
simulations have been done with to„ = nid and electric 
charge turned off. From these simulations one extracts 
physical results by extrapolation in the up and down 
quark masses. The calculated hadron masses or de- 
cay constants are expressed as an expansion in quark 
masses using expressions derived from chiral pertur- 
bation theory. Once the coefficients of these chiral 
expansions, which are related to the low energy con- 
stants in the chiral Lagrangian, are determined then 
one has an overcomplete set of relations (because the 
number of hadronic observables are much larger than 
the input parameters) between experimentally mea- 
surable quantities and quark masses. Using these re- 
lations, extrapolations to the physical values of hadron 
masses or decay constants specify the physical quark 
masses. A self-consistent determination of the quark 
masses in terms of hadron masses or vice versa would 
validate QCD. 

The success of this program requires that three con- 
ditions be met. First, the simulations should be done 
with three flavors of dynamical quarks and the input 
masses for all three quarks (both in the update and in 
the construction of external quark propagators, i.e., 



sea and valence quarks) should be light enough to lie 
within the range of validity of the 0{p*) chiral La- 
grangian [^ElEj. This condition guarantees that the 
extracted chiral coeflicients are the same as in QCD 
and reliable. Second, the simulations should be done 
at lattice scales small enough that discretization er- 
rors can be neglected or can be removed by a reliable 
extrapolation of the data to a = 0. Third, simulations 
should be done on large enough lattices so that flnite 
volume effects are negligible. 

There are two ways in which the renormalized quark 
mass at scale ^ is defined using lattice simulations 
done at scale 1/a: 

m_R(^) Z„,{y.,a) m[a) 

Za (0 I dMt)J{0) I 0) 

^^^^^-^-^^^ - (0 1 p.(t) j(o) 1 0) • 

The first method is based on the vector Ward iden- 
tity and mia) is the bare lattice mass. The sec- 
ond method exploits the axial Ward identity and uses 
two-point correlation functions with source J having 
pseudoscalar quantum numbers. For lattice formula- 
tions with an exact chiral symmetry, e.g. staggered 
fermions, the two methods are identical. The connec- 
tion between results obtained in the lattice regulariza- 
tion scheme and some continuum scheme like MS used 
by phenomenologists is contained in the Z' s. Their 
calculation introduces an additional source of system- 
atic error in all quantities whose renormalization con- 
stants are different in the two schemes; spectral quan- 
tities (hadron masses) which do not get renormalized 
are an exception. I later discuss the quantitative effect 
on quark masses of the renormalization factors needed 
to connect lattice results to those in the MS scheme. 

In this talk I will analyze the state-of-the-art lattice 
data and discuss their reliability with respect to the 
following sources of systematic errors. 

• The number and masses of dynamical quarks 
used in the update of gauge configurations. 

• Chiral extrapolations to physical quark masses. 

• Continuum extrapolations to a = 0. 

• The uncertainty in the calculation of the renor- 
malization constant. In particular I will dis- 
cuss the difference in estimates of masses be- 
tween using 1-loop perturbative estimate, the 
one obtained in the RI/MOM scheme, and the 
fully non-perturbative one using the Schrodinger 
functional method. 
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3 State-of-the-art quenched results 

The state-of-the-art quenched results are summarized 
in Tabled Two recent results deserve some elabora- 
tion (I have not included results from domain wall [ 
[7] or overlap fermions as they are still preliminary 
and do not include a continuum extrapolation). 

(i) Results from the SPQcdR collaboration [ 9[ su- 
percede all previous estimates from the ROME group [ 
1101 which is why the latter are not included in 
Table ^ The new calculations improve on previous 
results in the following ways, (i) The SPQcdR cal- 
culations have been done using the non-perturbative 
0(a) improved Sheikholeslami-Wohlert action at four 
values of the coupling, /3 = 6.0, 6.2, 6.4 and 6.45. 
Over this range the lattice scale changes roughly by 
a factor of two (0.1 0.051 fermi), so a reliable ex- 
trapolation to the continuum limit has been carried 
out. (ii) The N^LO (4- loop) relation is used to con- 
nect the RI/MOM scheme to the MS scheme for the 
renormalization constants, as well as in the running of 
the masses to the final scale 2 GeV. (iii) Results using 
both the vector and axial Ward identity method have 
been computed and compared. 

(ii) The QCDSF collaboration is in the process of 
updating their 1999 estimates of the strange quark 
mass [El- They use the same methodology and the 
same values of coupling, f3 = 6.0, 6.2, and 6.4, as the 
SPQcdR collaboration [IH], so their results will provide 
a detailed consistency check. Their latest unpublished 
quenched estimate is r^mf^^ = 0.341(2) which trans- 
lates to mf^ « 100 MeV. 

Discussion: The estimates in Table ^ show a wide 
range of values that seem to depend on the lattice ac- 
tion, the quantity used to set the scale and the quark 
masses, and the renormalization constant. What I 
would like to emphasize is that this spread does not 
imply that the lattice calculations are in conflict, sim- 
ply that in the quenched approximation one does not 
expect consistent results and the spread is a manifes- 
tation of that. 

Within the lattice community it has been known for 
some time that, in the quenched approximation, there 
is a roughly 10% variation in estimates of the quark 
masses depending on the quantity used to set the lat- 
tice scale a (ro or Mp or Jk or Mm etc.). This was 
quantitatively demonstrated by Wittig at LATTICE 
2002 [ Ej who converted the results from the four 
best simulations by SPQcdR, JLQCD, CP-PACS, and 
ALPHA-UKQCD collaborations to a common scale 
set by rg. The result, shown in Table [21 is that these 
four estimates of m^, when extracted using a com- 
mon scale setting quantity, i.e., evaluating mg/ro in 
the continuum limit, show a much smaller variation 



and agree within errors. This indicates that the de- 
pendence on the fermion action, fitting procedures, 
statistics, and renormalization constants (perturbative 
versus non-perturbative) used in the calculations are 
much smaller effects. 

The other lessons we have learned from quenched sim- 
ulations are: 

• Quenched simulations do not give consistent es- 
timates of quark masses. Estimates depend on 
the hadronic states used to set the quark masses. 
For example set using Mk differs by 15 — 20% 
from that set using Mk* or M^ when the scale 
is set by Mp. Thus, quark masses are sensitive 
probes of the effects of dynamical quarks. 

• With non-perturbative results for Z's in hand we 
can evaluate how well 1-loop tadpole improved 
perturbation theory works to convert lattice re- 
sults to MS scheme. We find that 1-loop esti- 
mates work to within 5% for VWI method, i.e. 
Z,n for Wilson like fermions, and at about 10% 
for the AWI method . Given that the rest of the 
errors, once a common scale setting quantity is 
used, are of this order, the collapse of results in 
Table (21 to a roughly common value is not sur- 
prising. 

The bottom line is that quenched simulations have 
allowed us to refine the numerical methods, and to 
understand and quantify all other sources of errors to 
within 5%. So removing this approximation becomes 
the next step in obtaining precise estimates. 

4 Discussion oi Nf = 2 results 

The CP-PACS collaboration [ (23 El set the stage 
for large scale simulations with dynamical fermions 
by providing results that are of comparable quality to 
quenched simulations with respect to statistics, num- 
ber of quark masses used in the simulations, and in 
the number of values of lattice spacings used in the 
continuum extrapolations. Their estimates displayed 
a number of desired features. The most striking was 
that estimates of the strange quark mass from four dif- 
ferent methods using the axial and vector Ward 

identity definition of the quark mass and using either 

Mk or M^ to fix rus were in agreement. The 

numbers ranged from 86.9(2.3) to 90.3(4.9). 

The JLQCD collaboration [ il9., has provided another 
measurement of quark masses with two dynamical 
flavors that complements results by the CP-PACS 
collaboration [ 1201 ■ There are, however, a num- 
ber of technical differences in the two sets of calcu- 
lations. The CP-PACS calculation used the 1-loop 
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Action 
Renorm. 


m 

= {mu + TOd)/2 


ms{MK) 




scale 1/a 


JLQCD 


Staggered 


4.23(29) 


106(7) 


129(12) 


Mp 


(1999)[E] 


RI/MOM 










CPPACS 


Wilson 


4.57(18) 


116(3) 


144(6) 


Mp 


(1999)[[TS| 


1-loop TI 






132^4 




CP-PACS 


Iwasaki+SW 




lllt^ 


Mp 


(2000)['T^ 


1-loop TI 








ALPHA-UKQCD 


0(a) SW 




97(4) 




fK 


(1999) [ 15 


SF 










QCDSF 


0(a) SW 


4.4(2) 


105(4) 




ro 


(1999)[ITn| 


SF 










QCDSF 


Wilson 


3.8(6) 


87(15) 




ro 


(i999)[nni 


RI/MOM 










SPQcdR 


0(a) SW 


4.4(1)(4) 


106(2)(8) 




ro 


(2002) [ 9 


RI/MOM 











Table 1. State-of-the-art quenched results for quark masses. The labels for the action used are: 0{a) SW for non- 
perturbative 0{a) improved Sheikhholeslami-Wohlert (SW) fermion action, and Iwasaki for an improved gauge action. 
The renormalization factors are calculated using the 1-loop perturbation theory with Tadpole Improvement (TI) or the 
non-perturbative RI/MOM and Schrodinger functional (SF) schemes. The quantity used to set the lattice spacing a is 
listed in the last column with ro = 0.5 fermi obtained from the static force relation rdV (r) / dr\r=ro = 1-65 f I17|. All 
estimates are based on extrapolation to the continuum limit. ma{MK) and ms(M^) refer to the strange quark mass 
extracted using Mk or to fix it. 



Ref. 


m,(g') 


Q' 


F 


msiro) 


JLQCD 


106(7) 


Mp 


0.90(1) 


95(6) 


CP-PACS 




Mp 


0.86(2) 


98(2)(t^) 


SPQcdR 


106(2)(8) 


Mk- 


0.87(3) 


92(2)(7) 


ALPHA-UKQCD 


97(4) 


fK 


1.02(2) 


99(4) 



Table 2. Comparison by Wittig of m^MS, 2 GeV, Mk) from quenched simulations by the JLQCD, CP-PACS, SPQcdR, 
and ALPHA-UKQCD collaborations. Given are the original estimates and the scale setting quantity Q' , conversion factor 
F from Q' to ro, and the converted mass. 





Action 
Renorm 


m 


ms{MK) 


ms{M^) 


scale 1/a 
(GeV) 


JLQCD 


Wilson+SW 


3.22(4) 


84.5(1.1) 


96.4(2.2) 


Mp 


(2002)[[Tn| 


1-loop TI 








(2.22) 


CP-PACS 


Iwasaki-l-SW 


o 4c;+0.14 
■^■^0-0.20 


89±6 




Mp 


(2000) [EOI 


1-loop TI 






(a->0) 


QCDSF-UKQCD 


0(a) SW 


3.5(2) 


90±?o 




ro 


(2003) 


1-loop TI 








[1.9-2.2] 


QCDSF-UKQCD 


0(a) SW 




85(1) 




ro 


(2003) [02] 


RI-MOM 








[1.9-2.2] 



Table 3. Recent Nf = 2 results for quark masses. 
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Figure 1. Estimates of vris from the — 2 simulations 
by the CP-PACS and JLQCD (points on the finest lattice 
with a « 0.09 fermi) collaborations. Linear extrapolation 
to the continuum limit are shown for all four definitions of 
the quark mass. 



mean field improved value of csw i^i the fermion ac- 
tion, whereas the JLQCD uses the non-perturbative 
value. CP-PACS used the Iwasaki improved gauge 
action whereas JLQCD uses the unimproved Wilson 
(plaquette) action. CP-PACS had results at three 
values of the lattice scale (a « 0.22, 0.16, and 0.11 
fermi) while JLQCD provide data at a single point at 
a = 0.0887(11) fermi. (The fourth point in the CP- 
PACS calculation at a = 0.0865 fermi was used only 
as a consistency check because it has small statistics.) 

The JLQCD estimates spoil some of the nice consis- 
tency shown by the CPPACS analysis. In particular 
if one combines data from the two calculations, the 
extrapolations of AWI{Mk) and AWl{M^) estimates 
give « 88 MeV, whereas those from YWI{Mk) and 
VWI(M0) extrapolate to « 93 MeV as shown in Fig.d 
JLQCD quotes their AW1{Mk) value, 8A.5tl^^°, as 
their best estimate assuming this method has the 
smallest a dependence. The difference between the 
AWI{Mk) and YWI{Mk) values is taken as an esti- 
mate of the systematic uncertainty. I discuss these 
data further below. 

The QCDSF collaboration [||24j is in the process of up- 
dating their Nf = 2 estimate given in [|^. The piece 
of the calculation still missing is a non-perturbative 
evaluation of the renormalization constants. They 
should be finishing this calculation soon, meanwhile 
their unpublished estimate, using perturbative esti- 



mates of renormalization constants, is rn^^(2 GeV) = 
85(11) MeV. 

5 Continuum extrapolation 

I will use the data and fits in Figure ^ to illustrate the 
systematic uncertainty associated with the continuum 
extrapolation and the associated issues of renormal- 
ization constants and the partially quenched approxi- 
mation. Even though, as mentioned above, the point 
at a « 0.09 fermi is obtained with a different gauge 
and fermion action and therefore expected to have a 
different coefficient for the 0(a) errors, nevertheless, I 
have taken the liberty of making a common fit to qual- 
itatively illustrate how [in] sensitive the conclusions are 
to current errors in individual points and to a linear 
fit that extends all the way to a = 0.22 fermi. 

The extrapolated values are based on a linear fit to 
four points. Keeping a linear term is appropriate 
since 0(a) errors have not been fully removed from 
the action or the currents, however this does not mean 
that higher order corrections are unimportant. Look- 
ing at the fits it clear that more high precision data 
at smaller values of a are required to include/exclude 
higher order terms with a reasonable degree of confi- 
dence. Given the spread, JLQCD choose 84.5l|^^y*' as 
their best estimate of since AWI{Mk) values show 
very little a dependence. The different extrapolations 
are accommodated by associating a large positive sys- 
tematic uncertainty to the central value. 

It is interesting to note that the two estimates using 
AW I extrapolate to k, 88 MeV, whereas those using 
the VWI to K. 93 MeV. This suggests that the dif- 
ference is not due to using Mk versus but due 
to AW I versus VWI methods. There are two differ- 
ences between these methods that could account for 
this discrepancy. First, the renormalization constants 
are different for lattice actions that do not preserve 
chiral symmetry (one needs Za/Zp in the AW I and 
Zm = ^/Zs for the VWI); and second, there is an 
extra complication, in the case of the VWI method, 
coming from having to determine Kc, the critical value 
of the hopping parameter corresponding to zero quark 
mass. The problem of the determination of Kc using 
VWI from partially quenched simulations leads to an 
additive shift in estimates of quark masses. A compar- 
ative analysis of m„_d suggests that this issue leads to 
no more than 1 MeV uncertainty in estimate of quark 
masses, so I concentrate on the Z' s for the explana- 
tion. 

Repeating the CP-PACS and JLQCD analysis shows 
that most of the difference comes from the renormal- 
ization factor that connects lattice results to those in 
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the MS scheme at 2 GeV, i.e. Zywi ~ ^-^Zawi- An- 
other way of stating this is that the lattice values of 
TOs are roughly the same for the two methods; it is 
the connection between the two schemes which leads 
to majority of the difference. What we know from 
non-perturbative calculations of renormalization con- 
stants in the quenched approximation is that Zp is sig- 
nificantly (by about 10%) overestimated by tadpole- 
improved 1-loop perturbation theory, whereas Zs and 
Za are much better approximated. If we assume that 
the same is true in the Nf — 2 case, then correcting 
for this in the CP-PACS/JLQCD data would boost the 
AWI results by about 10% since ma — {ZA/Zp)m, 
and explain the difference between AWI and VWI 
results. In that case all four fits would extrapolate to 
ms(MS,2 GeV) « 93 MeV. In the absence of non- 
perturbative estimates for the renormalization con- 
stants, my conclusion, based on the CP-PACS and 
JLQCD results, is to take a flat distribution between 
84 — 93 MeV as the best estimate for rag . This range 
also incorporates the QCDSF-UKQCD estimates. 

Note that my reservation of ~ 10% uncertainty due 
to the 1-loop Z's does not apply to the four sets 
of quenched data analyzed by Wittig because none 
of those calculations use, simultaneously, the AWI 
method and 1-loop estimates for Z's. 

6 Nf = 3 results 

Simulations with three flavors of dynamical quarks, 
all with masses < m^, represent a qualitatively big 
step forward. The reason for this favorable situation 
is that the coefficients of the chiral Lagrangian deter- 
mined by fitting data for observables obtained at dif- 
ferent masses to the corresponding chiral expansions 
are the same as QCD |3 Thus, as long as the 
simulations are done within the region of validity of 
O(p^) xPT, we can extract physical results even from 
simulations done at 3 — 18 times md- 

Very recently preliminary results from simulations 
with three flavors were reported by the MILC [ I^HI 
and CP-PACS/JLQCD [ EHI collaborations at LAT- 
TICE 2003. These are summarized in Table H 

The CP-PACS/JLQCD collaboration [1201 employ the 
same analysis as in their Nf = 2 study [1201 and use 
an improved gauge action as well as an 0(a) improved 
Wilson quark action. Analysis of the AWI data give 
nis = 75.6(3.4), and the analysis of the VWI data 
is not complete as the determination of Kc is not yet 
under control. Their most accurate number is from 
the AWI(Af_R-) method and AWI(M^) gives a consis- 
tent value but with much larger errors. The difference 
between the two estimates is folded into the error es- 
timate. 



Two major issues remain with the CP-PACS/JLQCD 
results. These are (i) residual discretization errors as 
the calculation has been done at only one lattice scale 
and (ii) the use of 1-loop tadpole improved Z's. To 
address the first requires more data which is a matter 
of time. On the second issue my reservation, that 
the 1-loop perturbation theory overestimates Zp by 
~ 10%, resurfaces. If this reservation holds up then 
their estimate could be as high as 85 MeV. Based on 
these reservations, their numbers suggest the rather 
large range nis = 70 — 90 MeV. 

The MILC Collaboration results are obtained us- 
ing improved staggered fermions. (An older estimate, 
based on an independent analysis of a sub-set of this 
MILC data, was reported by Hein [W, at LATTICE 
2002.) I have listed the results under AWI even though 
for staggered like fermions (lattice fermions with a chi- 
ral symmetry) the AWI and VWI methods are iden- 
tical. A major step forward in the MILC analysis is 
they perform a combined fit to data for A/^ at two 
sets of lattices (coarse and fine) using a staggered xPT 
expression that includes discretization and taste sym- 
metry violating corrections 

(M^-;'°°f)2 

/i {rrix + rriy) 
16/i 

+ -j^[2Ls - L^){mx + my) 

+ '^{2Lq - Li){2mq + TUs) + a^C , 

where nix + is the sum of the masses of the two va- 
lence quarks. Even though fit to a complicated SxPT 
expression with 44-46 parameter, having very precise 
data allows them to extract the central values and the 
associated errors estimates reliably. Their best esti- 
mates are rriu^d ~ 2.7(6) and = 70(15) MeV. A 
very large part of the error comes from the fact that 
to the 1-loop estimate for Zm they assign an overall 
~ 20% uncertainty due to the neglected O(a^) terms. 
In my opinion this is a conservative estimate of the 
0{al) uncertainty, especially since the 1-loop coeffi- 
cient for the improved (AsqTad) staggered fermions is 
small, (<0.12as) [EHl- The authors are clearly keep- 
ing in mind the lesson learned from quenched unim- 
proved staggered fermions where that 1-loop pertur- 
bation theory underestimated Zm (and thus nis) by 
almost 30%. 

A concern with the MILC simulation is the lack of a 
"proof" that the staggered fcrmion action describes 
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Action 
Rcnorm 


m 


ms{MK) 
AWI 


scale 1/a 
(GeV) 


JLQCD 
(2003) [ESI 
MILC 

(2003) [EH 
MILC 
(2003) [EZI 


Iwasaki+SW 
1-loop TI 
AsqTad 
1-loop TI 
AsqTad 
1-loop TI 


2.89(6) 
2.5(1) 
2.6(1) 


75.6(3.4) 
66(1) 
68(1) 


Mp 
(2.05(5)) 

Mp 
(1.6) 

Mp 
(2.2) 



Table 4. Recent Nj = 3 results for quark masses. AsqTad is a perturbatively improved version of Staggered fermions 
which reduces "taste" symmetry breaking. The quoted errors in the MILC results include both statistical and those due 
to varying g* in the 1-loop matching between (1/a ^ 2/a). The entries in the last column give the quantity used to set 
the scale and its value in GeV. 



four degenerate flavors in the continuum limit. Fur- 
thermore, there is the potential problem of loss of lo- 
cality of the action when taking the square root and 
the fourth root of the staggered determinant to simu- 
late two plus one dynamical flavors. These issues are 
being investigated [ 29 now that all other sources of 
errors are understood, and the community is moving 
towards providing precision results. Unfortunately, as 
of now there is no airtight argument that settles these 
issues. 

Based on these two preliminary calculations, and if 
forced to quote a single number, my choice is nis = 
75(15) MeV. This estimate is certainly very exciting 
and provocative. Furthermore, with simulations at 
more values of the lattice spacing and with different 
fermion formulations coming on line, this exciting re- 
sult will soon be refined. 

As mentioned before, the power of Nf = 3 analysis, 
provided all quark masses are small such that 1-loop 
xPT applies, is that the chiral parameters are those 
of physical QCD. Thus, in addition to estimates for 
quark masses the MILC collaboration [EHl extract the 
Gasser-Leutwyler constants from their fit. In particu- 
lar they find that 

2L8-L5--0.1(l)(tJ) X 10-3. (3) 

This is significantly outside the range 

- 3.4 X IQ-^ <2Ls-L5< - 1.8 x 10"^ (4) 

acceptable for m„ = 0. The same conclusion has been 
reached by the OSU group ["SIT. In short, lattice re- 
sults do not favor the possibility that to„ = is the 
solution of the strong CP problem. 

7 rus from QCD Sum Rules 

Three types of sum-rules have commonly been em- 
ployed to determine light quark masses. They are (i) 



Borel (Laplace) transformed sum rules (BSR's); (ii) 
finite energy sum rules (FESR's); and (iii) Hadronic r 
decay sum rules (these r-decay SR are a special case 
of FESR) . 

The starting point for the pseudoscalar and scalar 
QCD sum rules arc the axial and vector Ward identi- 
ties 

d'^V;;" = (m, ~mu)i:su: 

and the corresponding integrated 2-point correlation 
functions, e.g. 

= I J d'^xe">-^{0\T{d''Al{x),d''A40)}\0), 

= {md + m^)^ij d''xe'''-^0\T{P\x),P{0)}\0). 

^(g^) and ^'5(9^) are analytic on the complex plane 
with poles and cuts along the positive real axis. They 
can be calculated using OPE and perturbation the- 
ory for large (say > sq) and away from the 
cut. They are related, through dispersion relations, to 
spectral functions, for example, p5{s) — Im^'5(s)/7r. 

Finite energy sum rules are based on the observation 
that the spectral function has singularities (poles and 
cuts) only along the real axis as shown in Fig. El The 
contour integral shown in Fig. El is zero so 

/ w{s)phadromc ds ^ (f) w{s)\lpQCD ds 

Jo 27ri J\s\=so 

The left hand side (discontinuity along the real axis) 
is evaluated using a combination of experimental in- 
put and modeling for the spectral function. The right 
hand side is evaluated using the OPE and perturba- 
tion theory for the dominant mass-dimension D = 
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Figure 2. The contour integral used in finite energy sum 
rules. The poles and cuts are along the real axis. Duality 
refers to the matching, on average, between the perturba- 
tive and hadronic ansatze. 



term. The scheme and scale (sq) used for the per- 
turbation defines the scheme and scale in which the 
quark mass is defined, ^/{s) are conveniently chosen 
weights designed to improve convergence. Since the 
OPE is expected to break down near the real axis at 
So, recent analyses have employed "pinched" weights 
like w{s) = (1 — s/sq)" that have a zero at s = sq. 
The resulting sum rules are called pinched FESR. 

In the Borel transformed sum rule one integrates the 
spectral function of the vector current along the real 
axis 

B[Il]oPE = / e^"^^'^ Phadronicds 

Jo 

/"OO 
J So 

For the strange scalar channel the l.h.s. is proportional 
to {nis — ■niuf'- In the OPE, the dominant D = {) term 
is evaluated using perturbation theory and defines the 
scheme and scale at which the mass is evaluated. The 
breakup of the integral on the r.h.s depends on a suit- 
able choice of sq- It has to be large enough that the 
integral of the perturbative estimate of the OPE (sec- 
ond term) is reliable and yet small enough that there is 
experimental data on the spectral function up to that 
point. Otherwise there is a large gap in which p can, 
at best, be modeled. Second, the answer should be 
independent of the Borel Mass M . One cannot choose 
M too small as the transform gives more weight to the 
low s/M"^ region. This is good on the r.h.s. but un- 
fortunately on the l.h.s. it enhances the uncertainty of 
the higher dimensional operators. At the same time 
we want Sq/M"^ > 1 so that the unknown contribution 
of the "continuum" to phadronic is suppressed. 



There are three important questions central to the re- 
liability of all sum rules analyses. 

• How well does the operator product expan- 
sion converge? Furthermore, are the non- 
perturbative corrections, like quark and gluon 
condensates, instanton effects, and neglected 
higher order terms in the OPE small? 

• How well is the perturbative expansion for the 
leading terms in the OPE known and how well 
does it converge at the scale sq? 

• How well is the hadronic spectral function de- 
termined through a combination of experimental 
data and modeling? 

With respect to these points two major improvements 
have occurred over the last five years. These include 

• The perturbative series for the scalar and the 
pseudo-scalar sum-rules are now known up to 
ttg (four loops) [EI]. 

• Better models of the hadronic spectral function 
have been developed that satisfy a number of 
consistency checks. 

A number of hurdles, mainly in our ability to deter- 
mine the phenomenological spectral function, remain. 

• In the pseudoscalar sum rule for rriu + rrid, 
the hadronic spectral function includes the 
masses and widths of the kaon, if(1460), and 
if (1830) resonances (to extract m„ -I- nid from 
the pion channel the corresponding states are 
7r,7r(1300),7r(1770)). What are not known 
are the decay constants of the X(1460), and 
i4r(1830) and their relative phase. Also, the Ktttt 
continuum is modeled using resonant forms with 
or without chiral perturbation theory modifica- 
tions. The prospects of new data to improve the 
spectral function are small. 

• In the scalar sum rule for ttIs , the spectral func- 
tion starts at the Ktt threshold. Also known 
are the masses and widths of the Kq{14:30) 
and K^{19b0) resonances. Below the K^{U30) 
threshold, the spectral function is fairly well de- 
termined using the K^s data and the Kt: scatter- 
ing phases. What is not known is whether there 
is significant phase variation near and above the 
ifg(1950) resonance. Prospects of improving 
p(s) from B-factories are marginal. 

Without going into details which will be presented in [ 
m, my conclusion, based on the Borel and finite energy 
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sum rules, is that the most complete analysis for 
from the scalar channel gives = 99(16) MeV [I32|. 
For the pseudo-scalar channel it is = 100(12) [\'6'6\. 

The T-decay sum rules utilizes data for the ratio of 
semi-hadronic to leptonic decay rate 

j^v/A^tj ^ r[T- ^ ^^hadronsyz/^^.j (7)] 

where V/A, ij denotes the flavor {ud or us) of the Vec- 
tor (V) or Axial (A) current. Experimental data gives 
access to the u, d and it, s spectral functions. The per- 
turbative series for the 1-1-0 part of the D = 2 r-decay 
sum rule are known only up to (three loops) . The 
part of the r-decay series is known to . Here and 
1 refer to the angular momentum of the hadronic part. 
A summary of current estimates of rus from hadronic 
r-decay sum rules is given in Table 151 [ I34| . 

The interesting feature to note, notwithstanding the 
many improvements, are that the results in Table [S] 
have been fairly constant over the last three years. The 
central value of has stayed in the range 115 — 120 
MeV if unitarity of CKM matrix is imposed, and 
rris = 100 — 105 MeV if the Particle Data Group values 
for I Vus I are used. Both estimates have errors of about 
20 MeV. In short, the results are very sensitive to the 
value of I Vus I, and the difference between the two esti- 
mates is of the same size as all the other uncertainties 
combined. 

Drawbacks of the r-decay sum rule are: (i) The 
Cabibbo suppressed hadronic r-decay data has not 
been separately resolved into J = and J = 1 con- 
tributions and (ii) there is large uncertainty in the 
perturbative behavior of the scalar component. Some 
progress in reducing a sub-set of these uncertainties 
was recently reported by the GPJSP collaboration [ 
I39| where they used a phenomenological parameteri- 
zation for the scalar and pseudoscalar spectral func- 
tion in the OPE. These new results, shown in TableEl 
are, nevertheless, consistent with previous estimates. 

In terms of future prospects, we expect significant im- 
provement in the measured r-decay spectral function, 
especially above the K*. Meanwhile, it is clear that, 
at least as far as the central value of is concerned, 
pushing it significantly below 100 MeV is disfavored 
by the sum-rules analyses. 

8 Lower bounds on quark masses 

Even though the sum rule and Lattice QCD esti- 
mates overlap within combined uncertainties, the cur- 
rent best estimate of the lattice result (75 ± 15) MeV 
is tantalizingly small. The following question is often 



raised. Does this lattice number violate rigorous lower 
bounds predicted from a sum-rule analysis? My an- 
swer is NO and I give a brief justification for this [ 

m 

The most stringent bound predicted is the "quadratic" 
bound obtained by Lellouch, de Rafael and Taron [ 40 . 
It predicts, assuming perturbation theory becomes re- 
hable by g = 2 GeV, that ms(MS,2 GeV) > 100 
MeV. The Achilles' heel of this analysis is that the 
perturbative expression that enters into the quadratic 
bound has very large coefficients 

25 

3J^o^2 - 2(J^i)2 = 1 + — a(Q2) + 61.79a2(g2) 

O 

+ 517A5a^{Q^) + ... 
= 1 + 0.83 + 0.61 + 0. 51 + ... .(5) 

where the second expression has been evaluated at 
Q = 2 GeV with a = ajn « 0.1. Pushing Q > 2.5 
GeV already lowers the bound to (MS, 2 GeV) > 80 
MeV and going to a still safer value of Q = 3 GeV 
where as/n = 0.086 gives ms(MS, 2 GeV) > 60 MeV. 
With such a poorly behaved series it becomes an arti- 
cle of faith as to what Q is considered safe, and there- 
fore what value to take as the lower bound. 

There are two other bounds for which the perturbation 
theory is reasonable even for Q > 1.4 GeV. These are 
the So > given in [EHI, and the "ratio" bound given 
in [^2- Assuming that perturbation theory is reliable 
at Q = 1.4, both these bounds provide a much lower 
value, i.e., m,{MS,2 GeV) > 80 MeV. 

My bottom line on these bounds is that we should be 
concerned only if lattice results go significantly below 
rris = 70MeV. So, the interesting question is which 
estimate will change over time? Will the sum rule 
estimate come down to match the lattice value of « 75 
MeV, or will the lattice result rise to match the sum- 
rule value ms>100 MeV? Or will both change within 
their errors and come together in the middle? And 
finally, it will be interesting to understand which, if 
any, systematic error is being underestimated in the 
two methods. 
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106 ± 32 
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